Imaging the 3-D cosmological mass distribution with weak gravitational lensing 
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I show how weak gravitational lensing can be used to image the 3-D mass distribution in the 
Universe. An inverse relation to the lensing equation, relating the lensing potential evaluated at 
each source to the full 3-D Newtonian potential, is derived. I consider the normal modes of the 
lensing problem and clarify the equations using a small-angle approximation. Finally I consider 
the prospects of using this method to estimate the 3-D matter distribution from a realistic galaxy 
lensing survey. 
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I. INTRODUCTION 

Gravitational lensing provides us with the most direct 
and cleanest of methods for probing the distribution of 
matter in the universe 11,11. The lensing effect arises 
from the deflection of light by perturbations in the met- 
ric. These deflections stretch and contract bundles of 
light rays, causing the distortion of background galaxy 
images. Hence gravitational lensing does not depend 
on any assumptions about the state of matter. These 
distortions manifest themselves as a shear distortion of 
the source galaxy image , or a change in the surface 
number density of source galaxies due to magnification 
[^-H and can be used to map the two dimensional pro- 
jected matter distribution of cosmological structure. As 
the matter content of the universe is dominated by non- 
baryonic and non- luminous matter, gravitational lensing 
is the most accurate method for probing the distribution 
of this Dark Matter. Imaging of the Dark Matter distri- 
bution is a vital key to understanding its nature. 

The measurement of the gravitational mass distribu- 
tion in clusters of galaxies using gravitational lensing 
shear and magnification effects is now a well established 
technique, while on larger scales the detection of a cosmic 
shear signal ||] shows that the cosmological matter distri- 
bution can also be probed this way. A problem with these 
methods is the limited use of depth information, such as 
spectroscopic or photometric redshifts. Usually the inter- 
vening lensing matter distribution is approximated by a 
sheet. Some depth information can be introduced by lens 
tomography , where the background source galax- 

ies are divided into bins and the matter distribution can 
be approximated by a series of sheets. What is lacking is 
a method for estimating the full 3-D matter distribution 
from gravitational lensing. Here I address this problem. 

II. METHOD 

The metric of a perturbed Friedmann-Lemaitre- 
Robertson- Walker universe in the conformal Newtonian, 



or longitudinal, gauge is 

ds^ = -(1 + 2^)dt^ + a^(i)(l - 2^)dndr^ (1) 

where $ is the Newtonian potential, a is the cosmological 
scale factor, and we have assumed a spatially flat universe 
for simplicity. The Newtonian potential is related to the 
matter density field by Poisson's equation 

= A^GpraSa" - ^X-^^^ma-H, (2) 

where = 1/i/o ~ SOOO/i-^Mpc is the Hubble length, 
and r2„j is the present-day mass-density parameter. 

The lensing potential, <j), for a source in a spatially flat 
universe at distance r is given by 

4>{r)^2 j\r' (^—f^<^{r'), (3) 

measured at an angular position r on the sky. We have 
assumed the Born approximation, where the light path 
is unperturbed. This equation shows that the lensing 
potential is a radial projection of the 3-D gravitational 
potential, with a radial Greens function G(r, r') = 2(r — 
r')/rr'Q{r — r'), where Q{r) is the Heaviside function. 
The inverse relation to equation is 

<^{r)^^dydr^{r) (4) 

where dr — r.V is the radial derivative, and we assume 
the lensing potential has been appropriately smoothed to 
allow differentiation. This can be verified by substitution 
into equation and integrating by parts. 

The lensing potential is not an observable. The observ- 
ables are the dimensionless, symmetric, tracefree shear 
matrix, 7^ , which describes the distortion of the lensed 
image and the magnification, fi, which describes the 
change in area. The shear metric is 

= {d,d, - ^S^^d'^ 0, (5) 
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where di = r{5ij — = ''(V^ — fidr) is a dimension- 

less, transverse differential operator, and = did'' is the 
transverse Laplacian. The indices — (1,2). In this 
expression we have assumed a flat sky. The magnification 
is given by ^ 



(1-k)2-72|-i w1 + 2k. 



(6) 



The second approximation holds for weak magnification, 
where the lens convergence, k, is defined by the 2-D Pois- 
son's equation. 



(7) 



In principle there is another second-rank tensor which 
can be formed from a scalar potential; 



(8) 



where is the 2-D antisymmetric Levi-Civita tensor, 
the brackets indicate symmetrization of the indices and 
4>B is a pseudo-scalar potential. The symmetric tensor 
B, and pseudo-scalar have odd parity and therefore 
cannot correspond to the parity invariant matter density 
field. Hence we expect B = = (pB, and can use this to 
investigate noise and boundary effects in lensing. 

The lensing potential can be estimated from the shear 
field by the generalised Kaiser-Squires |H| relation; 



1'. 



(9) 



where is the inverse 2-D Laplacian operator. In prac- 
tice the shear field is only discretely sampled by galaxies 
so we must smooth the shear field to perform the differ- 
entiation. This also serves to make the uncertainty on 
the measured shear field finite, since each source has an 
unknown intrinsic ellipticity. This smoothing need only 
be perpendicular to the light path. There may also be 
intrinsic alignments of galaxies due to tidal effects dur- 
ing their formation [|T^,0|. However, these appear to be 
small at large distances [|14[]. 

The observable shear and convergence allow us to mea- 
sure the lensing potential up to an arbitrary function of 
r, the radial distance; 



(r) +-0(0, 



(10) 



where (j) is the measured gravitational potential and V'('') 
is a solution to the equations 
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7/> = 0. 



(11) 



This sheet-like gauge freedom arises because at each dis- 
tance the shear and convergence define the potential only 
up to an arbitrary constant. In principle this can ar- 
bitrarily change as a function of distance. As the re- 
construction of the Newtonian potential requires radial 
derivatives we have to smooth in the radial direction. 



This also transforms V' from an arbitrary radial func- 
tion to one which is smooth on the scale of the radial 
smoothing radius. The radial gauge freedom we see here 
is related to the so-called sheet-mass degeneracy which 
also arises as a constant of integration when deriving the 
convergence from the shear |l5| . However it is important 
to note that the radial freedom we see here also arises if 
the convergence is used to estimate the lensing potential. 

Since we expect the true lensing potential field to re- 
spect statistical isotropy and homogeneity, the observed 
0-ficld must obey the equation 



d_ 

dr 



(12) 



where the angled brackets (• • •) denote ensemble averag- 
ing, or tangential averages over light-paths at the same 
distance. Substituting back into equation (0) we find 



{d,4)) 



(13) 



is an unbiased estimate of the Newtonian gravitational 
potential. The derivation of equation ( p^ ) is the main 
result of this paper, and demonstrates that the full 3-D 
Newtonian potential, and hence the matter-density field, 
can be reconstructed from weak lensing observations. 

In practice %p will not be a major problem for large 
surveys, since the boundary conditions used for the in- 
version can be used to set the mean potential at a given 
radius to zero. This is fine if the survey is large enough 
and the mean potential is zero, but if the angular size of 
the survey is small the potential may not average to zero. 



A. Normal modes 

The lensing functions can conveniently be expanded 
in normal modes, this time taking the curvature of the 
sky into account. We define the 3-D spherical harmonic 
modes of a field by 



drr^ / dl7 0(r)j,(fcr)r/„(f]), (14) 



where — dt/a{t) is the causal horizon. For a 
spatially open universe with hyperbolic geometry the 
usual spherical Bessel functions can be generalised to the 
hyper-spherical Bessel functions, ji{x) Xg{nK,x). 

The harmonic moments of the potential are related to 
the convergence by a 2-D Poisson's equation, yielding 0] 



K£™(fc) = ~^e{i + i)(i)i„,{k). 



(15) 



Similarly the shear field can be decomposed into tensor 
spherical harmonics (e.g. pq|). 



"fvif')^\ - k'dk }^je,n{kMkr)Y^l,.^^^{r), (16) 
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where 



yE 

{£7n)ij 



and where gij and "; " are the metric and covariant 
derivative on the 2-sphere respectively. The 2ern{k) are 
then related to the lensing potential field by 



1 /(£_(_ 2)1 



(18) 



along the light-path and so shows lensing probes struc- 
tures along the light-path with wavenumber k\\ ~ 1/r, 
where r is the distance to the source. If we fix the sources 
at a single distance from the observer, the lensing inver- 
sion equation reduces to the 2-D Kaiser-Squires relation 



2X% 
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5(fc|l=r"\fci) (24) 



and so for large distances or small structure, we find that 
the lensing signal increases as a function of the source 
distance as k. ^ . 



Similar expressions exist for the non-gravitational parity- 
violating shear term. 

The normal modes of the Newtonian potential are then 
related to the lensing potential by 



$£m(fc) 



where 
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dk' -f \)bD{k' -k)- ai{k, k')](Pem{k') 

(19) 
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a,{k,k')^-j^ dr{k'r)\uikr)jdk'r). (20) 

The first term of equation ([l9| ) is simply minus the lens 
convergence, which reflects the angular distribution of 
the matter field. The second term contains all of the 
distance information and shows that the radial modes 
are correlated, as lensing accumulates with distance. 



B. The small-angle approximation 

For small-angles these equations can be simplified. It 
is useful to consider plane-wave solutions where the in- 
version equation becomes 



1 
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(21) 



and now fj, = k.r is a cosine angle. Combining this with 
the plane-wave solution for the convergence. 
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2 ^ ^ 
and Poisson's equation we find 



S = k^Xjr 

3 " 



(22) 



(23) 



This relation can be readily understood if we consider 
modes parallel and perpendicular to the light-path, k\\ = 

kfj, and k± = k{l — fJ.'^y^'^ respectively. While k± is 
related to the tangential scale of the structure being 
probed, the radial mode, k^ arises from the integration 



C. Measurement of the mass-density field 

These relations allow us to estimate the accuracy for 
reconstructing the 3-D density field from gravitational 
lensing. The covariance matrix of Fourier modes of the 
density field, S{k), estimated from lensing is 



Cs5{k)^P5{k) 



9 yn„,ii-ii^) 



(25) 



where we have assumed the underlying density field is 
statistically homogeneous and isotropic; 



{S{k)6*{k')) = (27r)3p(fc)(5i5(fc - k'), 



(26) 



and 5o(k) is the Dirac delta function. The second term 
arises due to Poisson sampling of the shear field by source 
galaxies, where Cims ~ 0.4 is the intrinsic dispersion of 
galaxy ellipticities, and n is the density of sources. 

The uncertainty on a measurement of the matter power 
spectrum from a finite galaxy lensing survey is 



27r 



APa(fc) ^ 

Ps{k) VFdhTfcl^ff' 



(27) 



where V^s = j d^r[P{k)/{P{k) + N{k))Y is the ef- 
fective volume of the survey, where N{k) is the Pois- 
son noise, and we sample the matter power spectrum 
in logarithmic intervals of dlnfc. Equation ( |27| ) shows 
the fractional uncertainty per modes divided by the 
square root of the number of effective independent modes 
that will fit into the survey volume. For a finite sur- 
vey the angular part o f the integra l should have lim- 
its l/{kR) < n < - l/ikRf and k > V^/R, 
where R is the size of the survey. For small scales 
this is dominated by shot-noise and fc^V^ff (fc)/27r^ « 
27/(56^)(r^P/e?^J2(fcAH)2(i^/Aff)lO(a„/a)^ which is 
strong function of the volume. 

Figure 1 shows the 3-D matter power spectrum and 
the expected uncertainty from an idealised galaxy lens- 
ing survey. The dotted line is the linear power spectrum, 
while the solid line is the nonlinear power JTzt . The model 
is a ACDM model, with J^a — 0.7 and — 0.3, nor- 
malised to the present-day cluster abundance pSl. The 
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FIG. 

A™(fc) = ^fc3p(fc)/27r- 



dimensionless matter power spectrum, 
for a ACDM cosmology, with 



57a ~ 0.7 and Qrn. ~ 0.3, normahsed to the observed abun- 
dance of clusters. The dotted line is linear theory, the solid 
line is nonlinear theory. The dot-dashed lines are the ex- 
pected uncertainty on a measurement of the power from a set 
of large-scale galaxy lensing surveys. 



upper dot-dashed line shows the uncertainty for a sur- 
vey with source density n = 10~^ /i~^Mpc~'^ and volume 
V = 10^ h-^Mpc^, similar to the Sloan Digital Sky Sur- 
vey. The lower lines are for larger surveys with the lower 
line for a survey ten time larger. Clearly the 3-D matter 
power spectrum can be recovered to good accuracy. 



III. SUMMARY 

I have derived an exact expression for the recon- 
struction of the full three-dimensional matter density 
field from measurements of weak gravitational lensing 
of source galaxies. This can be applied to matter fields 
in both the linear and nonlinear clustering regime, since 
the method is based on a new relation between the New- 
tonian potential and the gravitational lensing potential. 
The lensing potential can be estimated from observable 
lens shear or convergence fields. A new sheet-like de- 
generacy arises due to a constant of integration of the 
observable fields which is an arbitrary function of dis- 
tance. I have shown that this function can be removed 
by averaging over light paths. 

I have derived the normal modes of the lens fields and 
the reconstruction equation. Using a small-angle and 
plane- wave approximation I have derived a simplified set 
of equations. The small-angle approximation allows us to 
estimate the accuracy to which a realistic galaxy lensing 
survey can reconstruct the full matter distribution. I find 
that the reconstruction of the full distribution is possible 
with current surveys, if the data quality is sufficiently 
good. 3-D imaging is an interesting challenge for future 



lensing surveys. Since the addition of depth information 
adds an important new aspect to lensing, estimation of 
statistical quantities such as the 3-D matter power spec- 
trum should be possible with near-future surveys. Com- 
bined with the distribution of galaxies from the same 
survey, this will lead to direct constraints on the environ- 
mental aspects of galaxy formation, as well as a direct 
probe of the cosmological distribution of Dark Matter, 
which in turn will place strong constraints on theories of 
structure formation and the nature of the Dark Matter. 
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